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^\ ' Abstract. We present near-IR interferometric measurements of the K1.5 giant Arcturus (a Bootis), obtained at 

\^ , the IOTA interferometer with the FLUOR instrument, in four narrow filters with central wavelengths ranging 

' from 2.03 fim to 2.39 /jm. These observations were expected to allow us to quantify the wavelength dependence 

of the diameter of a typical K giant. They are compared to predictions from both plane-parallel and spherical 
model atmospheres. Unexpectedly, neither can explain the observed visibilities. We show how these data suggest 
the presence of a companion, in accordance with the Hipparcos data on this star, and discuss this solution with 
respect to Arcturus' single star status. 
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ryj , 1. Introduction excellent opportunity to study the deviations from a non- 

wavelength-dependent uniform-disk model (which is usu- 
Early K giants are often used as calibration sources for ally used to model the ca i ibrat i n sources, a step necessary 
photometry, spectroscopy and interferometry in the (near) m the ca i ibration process), and hence to investigate the 
IR because they offer a good compromise between bright- need for sophisticated models to calibrate interferometric 
g . ness at these wavelengths and compactness of the atmo- measur ements if high-accuracy visibilities are sought. Such 
sphere. The latter quality guarantees the absence of ex- a detailed study of Arcturus w as already performe d in the 
otic behaviour in both the line/band formation (deviations visible wave l en gth regime by iQuirrenbach dTafl il993) , 
from Local Thermodynamic Equilibrium (LTE), compli- who found good agreement with theoretical wavelength- 
cated temperature distributions, ...) and in spatial struc- dependent limb-darkening predictions. In this paper, we 
t ure - investigate the near-IR part of the spectrum. 

One such very popular star is Arcturus (K1.5-2III, 
a Bootis) . This star has been a stellar standard for We start out by discussing the interferometric calibra- 
many decades, not only in the infrared for spectro- tion process (Sect. 0). Then we present the data obtained 
sco pic /photometric wo rk, (e.g. as an ISO-SWS calibra- on Arcturus (Sect. |3J), after which we present different at- 
tor, iDecin et alJl2003al) but also as an IAU radial velocity mosphere models with which one can interpret the data 
standard fe.g. |Pearcelll955|) . A lthough its angula r diame- (Sect.0J. 
ter of about 20.20 ± 0.08 mas llPerrin et alj fl998i implies 

that it is too resolved by most interferometric configura- Surprisingly, some data points are not at all consistent 
tions to be used as a calibrator (Sect. [2}, it provides an with the proposed models, and seem to suggest a binary 

nature for Arcturus. We discuss this solution and con- 
* Michelson Postdoctoral Fellow front it with our current knowledge of this star in Sect. [SJ 

** Postdoctoral Fellow of the Fund for Scientific Research, Finally, in Sect.® we summarize our results and the open 
Flanders questions. 
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2. The interferometric calibration process: the 
need for well known calibrators 

Interferometric observations measure the spatial coher- 
ence of a given source between two or more apertures. 
Any optical defect on the line of sight such as atmospheric 
turbulence or in the instrument such as polarization mis- 
matches or dispersion will degrade spatial coherence. 

The classical solution to overcome this issue of deco- 
herence is to observe reference stars used as calibrator 
sources. The observed degree of coherence measured on 
the calibrator (also called fringe contrast) is compared to 
the value expected from prior knowledge of the source 
characteristics. This defines the interferometric efficiency 
(also called transfer function) of the instrument at the 
time when the calibrator was observed. In single-mode in- 
terferometers for which turbulent phase has been filtered 
out, the interferometric efficiency is relatively stable. Yet, 
in order to monitor any instrumental change to achieve 
the best accuracy, calibrator observations are interleaved 
with science target observations. 

In an ideal world, a calibrator should be point-like in 
order to yield an expected visibility of 100% to within 
an excellent approximation. However, if we receive non- 
zero flux from an object then it must have a finite 
angular diameter. Furthermore, the requirement of IR 
brightness and the more or less Planckian energy dis- 
tribution of stars, favours cool giants as IR interfero- 
metric calibrators. Consequently in practice, calibrators 
are always slightly resolved albeit with large visibilities. 
But an accurate visibility estimate requires an accurate 
diameter estimate. Some sources have been measured 
interferometrically or with the lunar occultation tech- 
nique. But for most calibrator sources an a priori es- 
timate is required. It has to be based on spectroscopy, 
photometry and modelling. Examp l es of such studies 
can b e fou nd in ICohen et all lll996t Il999|) . iBorde et alJ 
ll2002i). andlMerand et alJ l)2004|) for the near-infrared, or 
IVan Boekel et alJ l|2004f> for the mid-infrared. 

These indirect techniques prove to provide excellent 
wide-band diameter estimates with accuracies as good as 
a few percent, but higher spectral resolution and/or an 
expected visibility of the calibrator well below 100% will 
require a better understanding of the calibrator diameters. 
The observations presented here were made to empirically 
quantify the wavelength-dependence of the diameter of an 
early K giant and to test whether theoretical atmosphere 
models can be used to compute the wavelength-dependent 
diameter of other interferometric calibrators. 



3. The observations 

Before presenting the new data, we briefly discuss the in- 
strument, the calibrator star and the data reduction strat- 
egy 



3.1. The instrument 

The observations were performed in May 2002 at the 
IOTA (Infrared-Optical Telescope Array) interferome- 
ter located at the Smithsonian Instituti on's Whippl e 
Observatory on Mount Hopkins, Arizona l)Traubl ll998'). 
Several baselines of the IOTA have been used to sam- 
ple visibilities at different spatial frequencies. The data 
have been acquired with FLUOR (Fiber Linked Unit 
for Optical Recombinat i on) in the version described by 
ICoude du Foresto et alJ l)l998(l . FL UOR is the precursor 
of the now well-known VLTI/ VINCI l|Kervella et alJl2000l 
l2003j) . FLUOR has four outputs (two interferometric and 
two photometric) which were f o cused on a Nicmos 3 ar- 
ray developed bv iMillan-Cabetl l|2000l) and operated with 
frame rates ranging from 500 up to 2000 Hz. 

Observations were carried out in narrow bands with 
filters specially specified for molecular b ands and the con- 
tinuum region of cool stars in K (e.g. lDecinll200nl) . We 
will discuss in Sect. 01 which bands are present in the K 
band in the case of Arcturus and how they are expected 
to influence the interferometric observations. The narrow 
band filters transmissions are plotted in Fig. [5] They are 
named K203, K215, K222 and K239 where the three digits 
characterize the central wavelength: 2.03, 2.15, 2.22 and 
2.39 /xm respectively. The two continuum filters K215 and 
K222 sample the maximum transmission region of the K 
band. The K203 (H 2 bands) and K239 (H 2 and CO 
bands) sample the edges of the K band where stellar flux 
is attenuated by the poorer transmission of the Earths at- 
mosphere due to the absorption by water vapor. Any loss 
of coherence due to this water vapour is taken into ac- 
count by measuring the interferometric efficiency in each 
band separatly. 

3.2. The calibrator star: HR5512 

Observations of Arcturus have been bracketed by obser- 
vations of the calibrator star HR5512 (M5III) whose di- 
ameter is estimate d to 8.28 ± 0.41 mas as explained in 
IPerrin et alJ l(l998h . Using an M5 giant is of course rather 
risky in the sense that it probably also has a slightly ex- 
tended atmosphere with a wavelength dependent diame- 
ter. And since it is partially resolved, the interferomet- 
ric efficiency we used might be biased. Theoretically, one 
can expect an atmos pheric extension of abo ut 4-8 % for 
a late M giant (e.g. IWittkowski et"al] l2004h . At a spa- 
tial frequency of 25 arcsec -1 , this would induce a variation 
with wavelength of the calibrator's visibility below 0.5 %. 
Nevertheless, this effect must be kept in mind when inter- 
preting the Arcturus data. 

HR5512 d oes not appea r in the spect roscopic binary 
catalo gues of iBatten et all (|l989() and IPourbaix et alJ 
l|2004|) . Hipparcos does not list it as a visual binary ei- 
ther, tho ugh it does get the flag "suspected binary". In 
response, Mas on et al T lll999l) used speckle interferometry 
to possibly resolve the system, but they found no compan- 
ion. 
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HR5512 is known to be a sem i-regular (SR) vari- 
able star (e.g. lPercv fc FleminJl99^) for which Hipparcos 
found variations with an amplitude of 0. 1 mag and a pe- 
riod of only 6.3 days l|Koen fe Everll2002|) . However, more 
recently, samples of these Hipparcos short-period SR's 
have been the subject of other de dicated photometric sur- 
veys, such as those perfo rmed by iKoen et all l|2002() and 
iKerschbaum et all <j200l|) . which could not confirm the 
presence of any variab ility with a period below 35 days 
in any of their targets. IKerschbaum et al ] (|200ih suggest 
some instrumental artefact on the side of Hipparcos. We 
are therefore confident that the brightness of our calibra- 
tor HR5512 has not chan ged significantly du ring the ob- 
serving run. Furthermore. IKoen et al.| l|2002() find that in 
most SR's, the brightness variations are primarily due to 
changes in temperature, and not in diameter. This sug- 
gests that the diameter variations of HR5512 during its 
pulsational cycle are well below the 2 % limit of the con- 
stant temperature scenario and hence do not affect our 
science observations. 

3.3. Data reduction 

Fringe co ntrasts have been derived with t he procedure ex- 
plained in lCoude du Foresto et al.l l)l997j) . The bias in vis- 
ibility est imates due to photon noise has been removed 
following IPerrirJ l)2003al) . Turbulent corrugations of the 
wavefront are cleaned by the single-mode fibers of FLUOR 
except for the diff e rentia l piston between two apertures. 
IPerrin fc Ridewavl l|2004l have shown that the piston- 
induced bias is smaller than the 0.1% level for typical 
fringe frequencies of a few hundred Hz as used in FLUOR. 
The expected visibility of the calibrator is computed at 
the time it was observed. The interferometric efficiency 
is then interpolated at the time when Arcturus has been 
observed. Division of the fringe contrast of Arcturus by 
the interpolated interferometric efficiency provides the fi- 
nal visibility estimate. Correlations in fringe contrast and 
transfer function estimates are taken into account in the 
computation of error bars. The whole ca libration proce- 
dure has been published in IPerrirJ ll2003bl) . The accuracy 
of visibility estimates measured with the FLUOR setup is 
usually on the order of 1 % for most sources and can b e 
as good as 0.2-0.4% (|Perrinll2003bl IPerrin et al.H20044 . 
In other words, all known biases are much smaller than 
the effects reported in this paper and error bars are well 
estimated. 

Table fj] lists all measured visibilities in May 2002, 
grouped per filter. 

4. Interpretation of the data: comparison with 
theoretical atmosphere models 

Typical of present-day optical interferometric data is that 
they do not allow an inverse Fourier Transform (to ob- 
tain an image) but only offer information when they are 
compared to a model of the source intensity distribu- 
tion on the sky. This implies that one needs to know in 



Table 1. All measured squared visibilities of the 2002 run, 
grouped per filter. PA is the position angle (counted from 
North to East) of the interferometric baseline projected 
onto the sky. 



Date 


Filter 


spat. freq. 


PA 


V' 1 


g{V z ) 


(MJD) 




(arcsec -1 ) 


(deg) 






52427 


K203 


37.70 


128.44 


0.1656 


0.0054 


52429 


K203 


37.27 


122.93 


0.1885 


0.0050 


52429 


K203 


37.32 


123.35 


0.1901 


0.0050 


52430 


K203 


26.68 


122.37 


0.4030 


0.0117 


52430 


K203 


26.69 


122.63 


0.3961 


0.0100 


52431 


K203 


37.62 


131.03 


0.1688 


0.0048 


52427 


K215 


35.01 


121.57 


0.2347 


0.0052 


52429 


K215 


35.25 


123.75 


0.2167 


0.0055 


52430 


K215 


25.18 


122.97 


0.4533 


0.0071 


52431 


K215 


35.46 


131.41 


0.2288 


0.0054 


52427 


K222 


34.02 


122.02 


0.2581 


0.0055 


52429 


K222 


34.25 


124.28 


0.2568 


0.0057 


52429 


K222 


34.30 


124.74 


0.2535 


0.0056 


52430 


K222 


24.44 


123.38 


0.5336 


0.0108 


52431 


K222 


34.36 


131.70 


0.2522 


0.0058 


52425 


K239 


42.09 


93.42 


0.0867 


0.0054 


52427 


K239 


31.08 


116.61 


0.2937 


0.0083 


52429 


K239 


32.00 


125.74 


0.2861 


0.0111 


52430 


K239 


22.74 


124.06 


0.5719 


0.0099 


52431 


K239 


31.94 


132.04 


0.2886 


0.0082 



advance, through previous (non-interferometric) studies, 
which type of model to choose. For Arcturus, a first step 
could be the fitting of a uniform disk model for each indi- 
vidual filter which yields the apparent diameters at the dif- 
ferent wavelengths. These values can then be compared to 
theoretical predictions. While this strategy is satisfactory 
for the study of stars with very extended atmospheres for 
which accurate model atmospheres are not yet available 
and for which the observed change of size with wavelength 
is large, it is not sophisticated enough for Arcturus, for 
which we expect di ameter variations w ith wavelength of 
only a few percent l(Decin et al.ll2003ah . Hence, we would 
like to use a fully self-consistent theoretical model atmo- 
sphere and compute from this a wavelength-dependent, 
limb-darkened synthetic source intensity distribution on 
the sky. 

4.1. Arcturus atmosphere models and diameter 
determination 

In this section, we present our dedicated MARCS model 
atmospheres and the resulting spatial intensity profiles. 
These are then converted into synthetic visibilities and 
compared with the observations. 



4.1.1. The atmosphere models 

Arcturus was used as a primary calibrator for the calibra- 
tion of the ISO-SWS (Infrared Space Observatory Short 
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Wavelength Spectrometer) l|Decinl l2000) and w as for these 
purp oses modelled using the MARCS-code llPlez et all 
Il992l and references therein). This code is aimed at the 
modelling of atmospheres of cool (giant) stars, allowing 
both plane-parallel (PP) and spherical (SPH) geometries. 
Basic assumptions underlying the models are the follow- 
ing: LTE, hydrostatic equilibrium and conservation of en- 
ergy for radiative and convective flux. The radiative trans- 
fer equation is solv ed using an ALI (Approximate Lambda 
Iteration) method <|Nordlundlll984h with the OS (Opacity 
Sampling) technique as the statistical way of treating 
sp ectral lines. 

IScnmid-'BurJkl (|l975l) and lScholTl (|l985h found that the 
PP approximation was no longer valid in cool giants when 
a measure of the extension of the atmosphere, d, defined 
to be 



d = 



(1) 



with R the linear radius and tr ss the Rosseland optical 
depth, was larger than 0.05. 

Recent studies of the rol e of sph e ricity i n oxygen-rich 
cool s tars include tho s e of [Schol d l|l985|) . iBessell et alJ 
l|l989h and IPlez et all ljl9921) . Compared to PP models, 
the radiation field in a SPH model becomes diluted in the 
upper photosphcric layers. This usually leads to a decrease 
in temperature of the surface laye rs which can i nfluence 
significantly the molecular opacity l|Mihalasll978L p. 650). 

SPH models have to deal with the added complexity 
of an explicit radial dependence of all variables, plus a 
more complex equation of radiative transport, solved for 
different rays toward the observer: 



p dl v 1 p? dl v 
n v p dr n v pr dp 



S v 



(2) 



with p = cos 9 the cosine of the viewing angle, k v the 
monochromatic extinction coefficient per unit mass, p the 
density, l v the specific intensity and S v the source func- 
tion. Computationally, th is equation is sol ved by using the 
single-ray approximation l|Nordhindl 1 984t for a set of par- 
allel rays, intersecting the atmosphere at different values 
for the impact parameter. In computer time, the differ- 
ence between PP and SPH models is usually less than a 
factor of three although convergence is not guaranteed in 
the latter case. 

The high quality of the fit between the resulting 
synthetic spectrum and both the observed ISO-SWS 
spectrum and high -resolution FTS spectra presented in 
iDecin et all |2003b1. suggests that these hydrostatic mod- 
els offer a good representation of the true atmosphere of 
Arcturus. Since this comparison was limited to the ISO- 
SWS wavelength region (2.38— 45/im), we searche d the full 
Arcturus FTS spectral atlas l|Hinkle et alJll995|) in the 
FLUOR bandpasses for peculiar spectral features. Lines 
are sparse (mainly CN and some atomic lines) and well 
spaced up to the 12 CO 2-0 bandhead at 2.29 pm. A com- 
parison between the FTS spectra (summer and winter) 
and a spectrum synthesized from our atmosphere model 



4363.00 

1.0%- 
0.8 
0.6 



wavenumber [cm J 
4360.91 4358.82 4356.73 4354.64 



E 

o 0.4 



0.2 



0.0 





2.2930 2.2940 2.2950 2.2960 
Wavelength [/i.rn] 



2.2970 



Fig. 1. FTS spectra (summer and winter, in grey) com- 
pared to our synthesized spectrum (black) at the 12 CO 
2 — bandhead near 2.29 pm. The good match confirmes 
the quality of our model in the FLUOR wavelength region. 
The discrepancy in the Titanium line is due to a poorly 
known log gf value. 



at this bandhead is shown in Fig. ^ Clearly the match is 
excellent, confirming the quality of our model. 

From the quality of the fit between synthetic and ob- 
served spectra, we find little reason to prefer a SPH model 
over a PP one, suggesting that the atmosphere is quite 
compact (indeed, d — 0.02 for the model). Nevertheless, 
we computed both PP and SPH atmospher es for a grid 
around the stellar parameters determined bv lDecin et alJ 
(2003a) for Arcturus (listed in TableEJ, since the effects of 
sphericity may still be detectable in the intensity profile on 
the sky and hence also in the visibilities. From these mod- 
els, we derived a spatial intensity profile for each OS wave- 
length point in the K band. The full OS wavelength grid 
of our model contains about 150,000 wavelength points, 
with a resolution of R~20,000 at 2 ^m. These points are 
chosen in such a way as to accurately sample the total 
opacity as a function of wavelength and is based on ex- 
tensive atomic/molecular linelists and continuum opacity 
sources. The full K band spectrum of this model is shown 
in Fig. [21 together with the FLUOR filter profiles. 



Table 2. The stellar parameters from IDecin et alJ 
(2003b) and references therein, and the limits of the grid 
in which we searched for an optimal fit to the interfero- 
metric data (£ t is the microturbulent velocity). 



Parameter 


Value 


Lower limit 


Upper limit 


T eff (K) 


4320 ± 140 


4250 


4500 


logff (cm/s 2 ) 


1.50 ±0.15 


1.00 


2.00 


[Fe/H] 


-0.5 ±0.20 


-1.0 


0.0 


Ct (km/s) 


1.7 ±0.3 
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ntensity Profile 




Fig. 2. The transmission curves of the 4 narrow band 
filters used in the FLUOR instrument on IOTA together 
with a synthetic K-band spectrum of Arcturus. For cool 
late-type stars, the filter centered at 2.03 fim probes a H2O 
band in the spectrum (not present in Arcturus). For a 
K giant, the filters at 2.15 and 2.22 /mi both probe the 
continuum, and the 4th filter at 2.39 /im covers the CO 
first overtone band (and possibly also H 2 0). 

The resulting variety in intensity profiles (with wave- 
length) is largest for the model with stellar parameters 
T e g = 4250 K and logg = 1 .00 and is shown in Fig. [3] and 
Fig. 0] To improve display of the profile at the limb, we 
plot intensity as a function of fj, = (I — r 2 /i? 2 ^Q.y) 1 ' 2 , 
where r is the projected linear distance from the center of 
the stellar disk. The FLUOR-bandpass-integrated profiles 
are ovcrplotted in black. Even in the PP geometry there is 
a noticeable variety in profiles, but through the integration 
over the bandpasses we lose most of that information. In 
the SPH geometry, the variety in profiles is much stronger 
and even after integration differences do remain, mainly 
in the outer regions of the disk in the CO and H2O prob- 
ing filter 1 In this case, this is due only to CO, since the 
model predicts no lines of water vapor (the photospheric 
temperature is too high for this molecule to survive). 

Because of circular symmetry, these intensity pro- 
files can now be Hankel-tran sformed into visibility curves 
l|Hanburv Brown et alJll974|) 

\r 1 \ Jo I i< J ) J o ( nx <t> V 1 - M 2 ) M d/x 

V{x) = — j , (3) 

Jo I(P>) V d/i 

where x — B/X is the spatial frequency (arcsec -1 ) with B 
the projected baseline and A the wavelength, / the inten- 
sity profile, Jo the zeroth-order Bessel function of the first 

1 We remark that high spectral resolu tion interferometric 
observations, e.g. with VLTI/ AMBER dPetrov et alJ |5oOO) 
will resolve the individual spectral lines visible between /j, — 
0.2 and 0.0. The difference in visibility curve between such a 
spectral line and the continuum will be much stronger than 
that between the different FLUOR filters discussed here. 



Fig. 3. Intensity profiles for a plane parallel model at 
T c ff = 4250 K and logg = 1.0, normalized to the maxi- 
mum intensity in the K-band at the center of the stellar 
disk. Every gray curve was computed for an OS wave- 
length within the K band. Clearly there is a wide range 
of shapes over the stellar disk, but integration over the 
bandpasses reduces the differences. Because the model is 
semi-infinite for all viewing angles, /1 — corresponds to 
a singularity. 

ntensity Profile 




0.0 0.2 0.4 0.6 0.8 1.0 

Fig. 4. Intensity profiles for a spherical model at T c ff = 
4250 K and logg = 1.0, normalized to the maximum inten- 
sity in the K-band at the center in the stellar disk. Every 
gray curve was computed for an OS wavelength within 
the K band. We now see even more variety than in Fig.|3 
mainly toward the stellar limb. The 2.39 /im filter differs 
significantly from the others in this outer region 1 . 



kind, and 4> the angular diameter. Here, the angular diam- 
eter corresponds to the image size, i.e. the outermost point 
of the model (tr oss = 10 -7 ). It is the only free parameter 
of the model and is determined by fitting the theoretical 
curves to the observed data. The results of this fitting are 
presented in the next section. 
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4.1.2. Diameter definitions 

Comparison between this image size and the different 
types of diameters generally used in the literature is not 
trivial. In the case of a uniform disk (UD), the visibility 
is given by 

= 

nxq) 

with x and <j> as m Eq-El an d J\ the first-order Bessel 
function of the first kind. Fitting a UD model to the data, 
we find the set of angular diameters listed in Table |3] 
These values are not compatible with a single wavelength- 
independent UD, and in addition the fit quality appears 
very poor in most cases (column 3 of Table|3J) , thus demon- 
strating the need for atmosphere models when precisions 
of « 1 % are reached on angular diameters. 

Table 3. Uniform disk (UD) angular diameters derived by 
fitting Eq. 0] to the narrow-band data taken in respective 
filters or altogether. The third column gives the reduced 
chi-square of the fit (chi-square per degree of freedom), 
denoted Xi- 



filter 


4> (UD, mas) 


xi 


K203 


20.77 ±0.18 


9.20 


K215 


20.91 ±0.15 


8.26 


K222 


20.44 ±0.16 


0.40 


K239 


21.12 ±0.25 


2.01 


all 


20.72 ±0.16 


6.14 



Alternatively, one may be interested in the limb- 
darkened (LD) diameter. This diameter is the actual phys- 
ical diameter of the star, but is only well defined for stars 
with a compact atmosphere (i.e. in case of a high surface 
gravity). Since the outermost radial point of our atmo- 
sphere model (-R TRosa= io-0 is supposed to correspond to 
the physical boundary of the star, the LD diameter is ac- 
tually the image size (<fi) in the case of a SPH model (cf. 
Eq. 3 and the definition of /i). This is not correct for the 
PP models because these are semi-infinite for all viewing 
angles (with a singularity at fi = 0), therefore the value of 
the intensity profile at /i = is an extrapolation from the 
last calculated /i point. This last calculated /i point de- 
pends on wavelength, because different wavelengths come 
with different numbers of rays in the model. Figures OH and 
01 demonstrate this difference between PP and SPH mod- 
els. Consequently, the image size listed in Column 2 of 
Table 0] does not correspond to i? TRosa= i -7 and can not 
be converted into a tr oss = 1 diameter as is done for the 
SPH models below. 

A physically relevant and well defined diameter is the 
one corresponding to the tr ss = 1 layer. It can be de- 
rived from the image size (of a SPH model) if one knows 
the ratio of the outermost radial point in the model to 
the tross — 1 diameter. This ratio is similar to what is 



called the "extension of the atmosphere" , d, as defined in 
Eq.^ but for the models presented here, the outermost 
point corresponds to tr OS s = 10~ 7 . This TRoaa= i is also 
presented in Tabled 

Table 4. Angular diameter determinations with plane- 
parallel and spherical MARCS models. Only a very small 
improvement on the chi-square can be achieved by increas- 
ing the extension of the atmosphere of the model (SPH 
model with logg = 1.0), and such a low surface gravity 
cannot be reconciled with the spectral features. 



model 


PP, 4250K 


SPH, 4250K 


SPH, 4250K 




logs = 1.5 


logs' = 1-0 


logs = 1-5 


4> (mas) 


21.32 ±0.19 


22.37 ±0.19 


21.84 ±0.21 


0TROM=1 


NA 


21.16 ±0.18 


21.19 ±0.20 


d (Eq^) 

xi 


0.030 


0.057 


0.030 


5.63 


5.38 


5.51 



* OUt tO TRoss = 10 



Figures |S1 and El show a comparison of our theoretical 
visibility curves with the observed narrow-band visibili- 
ties. For the PP model (Fig. EJ, all four visibility curves 
coincide. The best-fit angular diameter for T e g = 4250 K, 
logg = 1.50 and [Fe/H] = -0.50, is 4> = 21.32 ± 0.19 mas. 
The chi-square per degree of freedom or reduced chi- 
square is then xi — 5.63. 

For the SPH model with the same stellar parameters 
(Fig. [BJ, the filter centered within the CO band around 
2.39 fim shows a slightly larger star, though that is barely 
visible in the plot. We obtain (f> — 21.84 ± 0.21 mas for 
T cff = 4250 K, logg = 1.50, [Fe/H] = -0.50. However, 
the improvement is small since xi = 5-51 for this model. 
Some improvement can be achieved by reducing the sur- 
face gravity of the model (increasing the extension of the 
atmosphere): this leads to a significantly different visibil- 
ity curve for the CO filter, but the improvement in terms 
of Xr i s very small (cf. Table^J, and such a low surface 
gravity c annot be reconcile d with the observed spectral 
features (|Decin et al.ll2003b(l . 

We remark that our newly derived diameters are 
slightly (but significantly) larger than the LD di- 
ameters determin ed in K broadband (20.91 ± 0.08, 
Perrin et all Il998h and in the optical (21.0 ± 0.2 mas, 
Ouirrenbach et all 1996). The bias introduced by the 
use of a single calibrator is taken into account in the 
error determination and should thus not be responsi- 
ble for this discrepancy. The diameter of Arcturus was 
also estimated with spectro-photometr i c techniques (e.g. 
Cohen et all Il999t iDecin et alJ l2003bl: IVan Boekel et alJ 
2004|), but uncertainties on those diameters are generally 
much larger and therefore compatible with both our re- 
sults and those discussed above. The uncertainties on the 
distance (11.26 ± 0.09 pc) derived from the parallax, the 
surface gravity and the mass, make it impossible to ob- 
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Fig. 5. A comparison between the FLUOR visibility mea- 
surements of Arcturus and a plane-parallel MARCS vis- 
ibility model (T cff = 4250 K, logg = 1.50; \\ = 5 - 63 for 
cf> = 21.32 ± 0.19 mas). The four theoretical curves co- 
incide, while there is a significant dispersion of the data 
points. 



Fig. 6. For the spherical model with slightly lower surface 
gravity (T cff = 4250 K, log 5 = 1.00; \l = 5-38 for tf> = 
22.37 ± 0.19mas), the visibility curve in the CO-probing 
filter no longer coincides with the others, corresponding 
to a slightly larger diameter, but it still cannot match the 
measured visibilities. 



tain an accurate estimate of the angular diameter which 
does not depend on observed flux levels and direct (inter- 
ferometric) measurements. T he recently derived em pirical 
surface brightness relation bv lKervella et all l)2004|) yields 
a diameter of 21.2±0.2 (assuming a 0.02 mag error on the 
K magnitude) which is well in agreement with our result. 

4.2. Discussion 

To our surprise, neither PP nor SPH models can explain 
the observed visibilities. Indeed, the probability to obtain 
« 5.6 with 19 degrees of freedom is as low as 10 -14 , 
and there appear to be systematics in the residuals. In the 
following, we investigate possible sources of these discrep- 
ancies. 

An important characteristic of the residuals is that 
they not only point at a problem with the wavelength de- 
pendence of the visibility, but also at a problem with the 
shape of the visibility curve which is not consistent with 
that of a limb-darkened disk: the large discrepancy at a 
low spatial frequency (25 arcsec -1 ) suggests a structure 
at least a few times larger than the stellar disk. 



4.2.1. Calibration problems 

A first check that needs to be done concerns the cal- 
ibrator source: as discussed in Sect. the extension 
of the calibrator's diameter should contribute at most 
only 0.5% to the residuals at 25 arcsec -1 . It is however 
possible that HR5512 shows a more complicated atmo- 
sphere/circumstellar environment, but some simple calcu- 
lus shows that ignoring this effect in the calibrator actually 
leads to an underestimation of the variation in visibility 
with wavelength in the science target: if the extension of 
the calibrator is affecting the calibration, then the prob- 
lem for Arcturus is actually even larger than the few per- 
cent reported here. Only if we assume the data points that 
fall below the fit (e.g. those at 26 arcsec -1 ) to be correct, 
can the extension of the calibrator cause the other points 
to be overestimated. However, these data points are not 
at all compatible with p hotometric d iameter estimates of 
Arcturus (e.g. 20.8 mas, |Decin 2000J). Moreover, in such 
a scenario, we expect problems mostly in the wavelength 
dependence and not in the shape of the visibility curve. 

Conclusive evidence against problems with the cali- 
brator is provided by a limited set of similar narrow- 
band FLUOR data on Arcturus from a run in May 2001 
(Table EJ). Unlike the 2002 data presented in this paper, 
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Spatial Frequency (arcsec ) 

Fig. 7. Narrow-band data on Arcturus obtained in 2001 
with different calibrators. The best fitting LD disk model 
(full line) shown here (xj? = 64.5) leaves the same kind of 
residuals as seen in the 2002 data. 

these 2001 data were calibrated with different calibrators, 
not including HR5512. They are presented in Fig. to- 
gether with the best fit LD disk model. It is clear that 
these data show very similar residuals. We thus conclude 
that the effects reported in this paper cannot be due to 
the calibrator. 



Table 5. All measured squared visibilities of the 2001 run, 
grouped per filter. 



Date 


Filter 


spat. freq. 


PA 


V' 1 


<t(V) 


(MJD) 




(arcsec - ) 


(deg) 






52043 


K222 


23.67 


98.30 


0.5079 


0.0054 


52047 


K222 


45.27 


95.51 


0.0684 


0.0021 


52043 


K215 


24.34 


95.47 


0.4522 


0.0044 


52045 


K215 


45.27 


89.02 


0.0689 


0.0020 


52047 


K215 


46.60 


93.07 


0.0525 


0.0022 



The largest deviations from the expected visibility 
where observed on only one night (MJD=52430), and one 
might wonder whether observing conditions showed any 
peculiarities that night. We found no sign of this night 
being any different from the others of that run. Strong 
piston often induces a broadening of the power spectrum 
of the affected fringe scans, but this is not seen in the data 
presented here. 

4.2.2. The photosphere of Arcturus 

Recently. iRvde et alJ (|2002l) reported the presence of H 2 
lines in a high-resolution Texes spectrum of Arcturus, 
which were not present i n the MARCS mode l used for com- 
parison (also used hereL lRvde et all l|2002|) show how the 
presence of these lines can be produced in the synthetic 
spectrum by imposing a slightly lower temperature in the 



outer layers (logr 50 o < —4) of the model-atmosphere. 
Although this new temperature distribution will influence 
somewhat the location of the line-forming region, the fact 
that the differences with the original MARCS model are 
only visible in a high-resolution spectrum suggest that 
it should not influence our low-resolution interferometric 
measurements. 

Still, recent results on supergiant stars demonstrate 
that additional, non-photospheric layers of water can have 
such a combination of geometry and temperature that the 
extra absorption of stellar flux is filled in by the emis- 
sion of the outer regions of the layer, leaving little trace of 
the significant additional water column density at medium 
spectroscopic resolution (e.g. lOhnakal l2004h . This ad- 
ditional water opacity might still be detectable through 
individual lines in high-resolution spectra, such as the 
Texes spectrum and will significantly influence the inter- 
ferometric observations. However, these additional water 
layers do create a pseudo continuum, affecting the spectro- 
photometry flux levels within the molecular bands. While 
this effect might still be minimal for one given band, the 
same temperature-geometry combination of the layer will 
result in a strong sig nature in another band. For Arcturus, 
iDecin et al.l l)2003bj) show that discrepancies between our 
MARCS model and the ISO spectrum are not present in 
any water band between 2.38 and 12 /im. Moreover, the 
existence of non-photospheric molecular layers is generally 
linked to the pulsation of the central star, and is therefore 
unlikely in the case of Arcturus. 

We might also wonder whether the pressure scale 
height coming out of the model is correct. It is known 
that this is not the case for cooler and more luminous 
stars llPerrin et alJl2004b(l . However, fitting a model with 
a larger atmospheric extension did not improve the \ 2 sig~ 
nificantly (see Tabled} , and the decrease in surface gravity 
required to obtain such an atmosphere is not compatible 
with spectroscopic data. 

Furthermore, we must note that the dispersion seen in 
the FLUOR data is not only a molecular-band vs. con- 
tinuum effect: the data obtained on the shortest baseline 
(around 25 arcsec -1 ) show that even for the continuum 
filters, not all points are consistent with a limb-darkened 
disk geometry. This feature actually rules out the hypoth- 
esis of an under-estimation of the extension of Arcturus' 
atmosphere: even if more pronounced than in the SPH 
models, it could never influence the intensity profile in 
the continuum to this extent. 

Another hypothesis could be the presence of spots 
on the stellar surface. However, this extra small scale 
structure would produce oscillations in the visibility curve 
whose period must be larger than the first null spatial fre- 
quency and hence it cannot explain our data. If we believe 
these data points are correct (and we see no reason not 
to), more exotic solutions should be explored. 

We are of course aware of Arcturus' status as cali- 
brator source for many different instruments/techniques 
and of the research done on this object in that framework 
but the following subsection will show that some further 
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investigation is nevertheless necessary to explain our in- 
terferometric data and, for example, also the Hipparcos 
data on Arcturus. 

5. The binary hypothesis 

In this section, we investigate the possibility that the dis- 
crepancy between the MARCS model and our interferomet- 
ric data might be due to the presence of a faint companion. 
Moreover, we discuss the consistency of this hypothesis in 
the light of previously published works. 



Table 6. Best- fit parameters with formal errors for the 
binary model. 





This work (K band) 


Hipparcos (V band) 


p 


212.7 ± 1.5 mas 


255 ± 39 mas 





157.6 ±1.7 deg 


198 deg 


r 


50.1 ±5.6 


21.5 ± 5.9 


Am 


4.25 ±0.12 


3.33 ±0.31 


<Px 


21.56 ± 0.05 mas 




02 


(fixed) 




X? 


2.6 





5.1. Binary model 

Let us denote by Vi and V2 the visibilities of Arcturus and 
its companion considered as a single stars. The companion 
is assumed to be unresolved {<p2 — 0), therefore V2 = 1, 
and for V\ we use the limb-darkened model of Eq. [3J The 
contrast ratio between the two stars, r = F1/F2 with F 
the received flux, is assumed to be identical for all four 
filters. This might be too crude an approximation if the 
spectral types of the two stars differ significantly, but we 
show in Sect. l5.2lH that this cannot be the case. We neglect 
the motion of the companion during the 7 days of our 
observation run, so that its separation p and position angle 
9 would remain constant. The squared modulus of the 
visibility for the binary then reads 



V 1 



% 2 



v 2 2 



2rViV 2 cos[2 7 r( J B/A) / 9Cos(6' - 



(r ± l) 2 



(5) 

where (B,9b) are the polar coordinates of the interfero- 
metric baseline vector projected onto the sky. Note that 
for this equation to be valid, the companion's fringe packet 
should be well overlapped with the primary star's fringe 
packet. Since the narrow filters used here guarantee about 
55 fringes in the fringe packet, while the derived separa- 
tion between the 2 components (Table corresponds to 
only 7 fringes, Eq. [SJcan be used. The parameters p, 9, r 
and (pi are estimated by minimizing a standard \ 2 with 
the Levenberg-Marquardt procedure. One should be par- 
ticularly careful in the search for the global minimum as 
trigonometric functions in Eq. [S] cause the \ 2 hypersur- 
face to have many local minima. To overcome this prob- 
lem, we use a four-dimensional grid of initial guesses to 
run the minimization, and keep the best x 2 of all. With 
this set of data, the estimation is made more difficult by 
the fact that the visibility curves were not sampled with 
the idea of looking for a binary (a continuous sampling at 
one wavelength would have been much more efficient). 

The best solution (Table EJ has a reduced \ 2 °f 2.6, 
twice as low as the best single-star model. For this solu- 
tion, we have plotted in Fig. the expected evolution of 
the squared visibility during our five nights on Arcturus 
together with the measured values. For comparison, Fig. [5] 
presents the same thing for the single-star SPH model. 
From the marginal x 2 curves in Fig. 1101 we see that there 
is no ambiguity on the global minimum with respect to r 
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Fig. 8. Squared visibility as a function of time, and the 
best fitting binary model. 
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Fig. 9. Squared visibility as a function of time, and the 
best fitting single-star SPH model. 



and <pi, the situation being more delicate for p and 9 ow- 
ing to the trigonometric functions in the model. Figures lTTl 
and 1121 show the confidence regions for the best solution 
as contours of the \ 2 surface in the two subspaces (p, 9) 
and (r, (pi ) . The form of Eq. [5] favours the correlation be- 
tween the parameters in each of these pairs (98 % and 
57% respectively), something that can be seen from the 
elongated shapes of the confidence regions as well. 
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Fig. 10. Marginal \ 2 f° r separation p, position angle 8, 
contrast ratio r and primary angular diameter (f>\ . Global 
minima are marked by a dotted line. 




Fig. 11. Confidence intervals as a function of separation 
p and position angle 0. 
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Fig. 12. Confidence intervals as a function of contrast 
ratio r and primary angular diameter <pi. 



5.2. Discussion 

Although Arcturus is generally believed to be a single 
star, a literature study reveals that Hipparcos also lists 
Arcturus as a do uble star This result caused some debate 
at that time fe.g. lGriffinlll998HTurner et al.lll999j) and in 
response, the Hipparcos team re-analyzed the data, and - 
though they arrived at th e same conclusion - flagged th e 
solution as "unreliable" l|Soderhielm fc Mignardl Il99sfl . 
Surprisingly, our solution agrees quite well with the 
Hipparcos solution (Table and we are therefore con- 
vinced that this hypothesis deserves further investigation. 
In the following, we will first rule out the possibility of 
a background star and a double calibrator, then we re- 
view briefly all relevant published works on Arcturus and 
finally confront these with our new results. 

5.2.1. Background star 

As a first step, we can rule out the possibil- 
ity of a visual binary, in the sense of a purely 
optical effect, since the high proper motion of 
Arcturus (—1093.43 a nd — 1999.43mas/yr in a resp. S, 
IPerrvman et alJ Il997|) would have increased the sepa- 
ration with the hypothetical background star by more 
than 10" since the time of the Hipparcos observa- 
tions, clearly separating the fringe packets of both 
stars. Furthermore, we obtained a new optical image 
of Arcturus and its neighborhood using the new CCD 
camera MEROPE on the 1.2 m Belgian Mercator tele- 
scope (http://www.mercator.iac.es/instruments). Images 
were made in all filte rs of the Geneva system (e.g. 
iRufener fc Nicoleilll988h but only in the (narrow) U band 
and with an exposure time below 0.01s did the flux not 
reach the saturation limit of the camera. This image con- 
firmed that no other star of sufficient magnitude is present 
at Arcturus' location in May 2002 (about 5" north-east 
from its current location) and at the time of the Hipparcos 
observations. 

5.2.2. Binary calibrator 

If the calibrator HR5512 were a close binary with a suf- 
ficiently bright companion, using a UD model to com- 
pute the interferometric efficiency might introduce the 
calibrator binary signature into the science observations. 
However, in Sect.|3|we already argued that HR5512 is not 
a known binary and in Sect. 14.21 we showed how fitting the 
2001 data set (for which other calibrator stars were used) 
with a LD model yielded similar residuals, resulting in 
a Xr = 64.5. We are therefore confident that the effects 
reported here are not due to binarity of the calibrator. 

Fitting a binary model to this 2001 dataset on 
Arcturus is possible, and returns parameters roughly con- 
sistent with those determined from the 2002 dataset: a 
separation of 188.5 ± 5.3 mas and a magnitude difference 
of 3.8 mag. The resulting x? of 1.65 is a clear improvement 
over the LD disk model. We must remark that the 2001 
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position angle (PA = 78°) differs significantly from that 
of 2002, but with such a limited dataset, this result is not 
very meaningful. 

5.2.3. The literature 

Radial velocities: The radial velocity of Arcturus has been 
studied for over a century now (e.g. lLordlll90^) and, be- 
cause of the stability of its velocity, ha s more recen tly 
become an IAU radial velocity standard l|Pearcelll955|h 

Nevertheless, velocity variations on several time-scales 
have been reported in the last decades. Ilrwin et alJ Jl989) 
obtained precise measurements of the radial velocity on 43 
occasions between 1981 and 1985. These show a range of 
500 m/s, with both short-period and long-period variabil- 
ity. For the long-period variability, they find an amplitude 
of 120 to 190 m/s and a period of 640-690 days. Since this 
period is longer than the fundamental radial mode of os- 
cillation, they discuss other possibilities: convection cells 
and dark spot are not really compatible with the absence 
of line- width changes. A companion with an Msini of 1.5 
to 7.0 Jupiter masses could explain their data. Explaining 
the long-period variability with beating of higher frequen- 
cies did not work well. In the last decade, a mode of pulsa- 
tion with an even longer period has been discovered: the 
gravity-mode. However, the periods associated with the 
g-mode are still well below 2 years, and it is not certain 
that these pulsations are observable through an extended 
giant atmosphere (Mazumdar, private communication). 

Although he does not presen t a det ailed analysis of his 
data in this context, ICochranl l|l988f) also notes a long 
period variability in his observations. Note that solar- 
like oscillations have recently bee n observed in Arcturus 
llMerlmelllflflrllRetter et a,l.ll2003h . 

Astrometry: No astrometric evidence for binarity was 
found by Hipparc o s dur ing its 4 year lifetime. However, 
iGontcharov et alJ <|200lf) combined astrometric ground- 
based catalogues containing epochs later than 1939 and 
the Hipparcos catalogue, to obtain new proper mo- 
tions and to detect non-linear astrometric behaviour. 
In their catalogue ("Proper Motions of Fundamental 
Stars"), Arcturus is listed as being an astrometric binary. 
Gontcharov (private communication) confirmed that the 
astrometric offsets (the residuals after subtraction of all 
linear motions) are significantly higher than what can be 
expected for a single star system. Periods of 5 and 20 years 
appear to be present in these data. 

Direct evidence: The only other direct evidence is the vi- 
sual detection of a companion by Hipparcos. The pub- 
lished results are listed in TableEJ Originally, no relative 
motion of the system was detected during the Hipparcos 
lifetime. However, after re-analysis of the data, a small 
relative motion of 4-8"/yr was detected, though at a PA 
which is "almost at right angl es" with the published value 
( Soderhielm & Mitrnar d1ll998l) . 



All other attempts at a direct detection of a compan- 
ion around Arcturus have returned no positive results. 
These include the non-detection in the Ha filter with the 
AO system on the Mount Wilson 100 Inch Telescope by 
iTurner et alJ l)l999l) . However, looking at Fig. 1 of this 
paper, a companion at a 255 mas separation would be 
located on the 2.0 x 10 3 contour of the primary's Point 
Spread Function (PSF). Assuming that the PSF peaks at 
3.0 x 10 4 at least, a contrast ratio of 20, i.e. 1.5 x 10 3 
at most, would be below the level of the PSF, making 
the detection difficult. Moreover, the reconstructed image 
is not diffraction-limited (as 65 mas would be the size of 
the first Airy ring), and not circularly-symmetric, showing 
that residual aberrations are present. Now speckle noise is 
certainly a concern for faint companion detection, so we 
believe these observations cannot definitely rule out the 
companion as found by Hipparcos. 

Nevertheless, another non-detection with Keck aper- 
ture masking, and the (apparent) absence of another fringe 
packet in IOTA/IONIC (broadband) H-band interfero- 
metric observations (Monnier, private communication) 
suggest Arcturus to be a single star indeed. Unfortunately, 
these data/results wer e not conclusive ( Monni er, pri- 
vate communication). iQuirrenbach et a"fl l|l996h found 
a good agreement between their optical interferometric 
data on Arcturus and theoretical limb-darkening profiles. 
Especially the very low visibility measured around the first 
null (their Fig. 3c) seems incompatible with an unresolved 
companion only 25 times fainter than Arcturus. 

Spectrophotometric observations It is noteworthy here 
that recent NLTE atmospheric modelling and spectrum 
synthesis calculations have been unable to match the ob- 
served spectrophotometric data in the violet and near UV 
bands: the flux observed between 300 and 400 nm amounts 
to on ly half of the predicted flux l|Short & Hauschildil 
l2003h . While this unresolved discrepancy might have noth- 
ing to do with binarity, it is a further clue that there may 
be more to Arcturus than is known. 

5.2.4. Nature of the companion 

According to our binary model (Table EJ), the contrast ra- 
tio primary/companion in the K band amounts to about 
50. The UV problem mentioned in Sect. I5.2~3l not taken 
into account, there is no clear photometric evidence of a 
composite SED/spectrum (Fig. 1131) . Therefore, the com- 
panion should emit as an object of fairly similar spectral 
type: if it were significantly warmer, the contrast ratio in 
the optical would be much lower, resulting in a clearly 
composite SED. If the companion were cooler (but co- 
eval), it would have to be a giant further in its evolution 
and thus more massive and brighter than the primary 2 . 
This is clearly not a possibility. We can therefore conclude 

Note that a Brown Dwarf is not a possibility for it would 
not reach the required brightness (4 mag) in the K band as 
derived from the contrast ratio r. 
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Fig. 13. Upper panel: Comparison of new UBVG 
(Geneva), JHKL (SAAO) and IRAS fluxes with a syn- 
thetic SED computed from our MARCS model atmosphere 
for Arcturus (A). There is no clear indication of a com- 
posite SED. Lower panel: the companion we suggest here 
(Kurucz model shown), would not be detectable in the 
residuals after subtraction of the primary. 



that this contrast ratio is a good upper limit on the lumi- 
nosity ratio of the two stars. In case of a slightly warmer 
companion, the luminosity ratio could possibly be a fac- 
tor of two lower, i.e. 25. Ado pting a total luminosity of 
196(±21)L Q as determined bv lDecinl |2000), and a lumi- 
nosity ratio range of 25-50, the primary would contribute 
188-192 L & and the companion 4-8 L©. 

A companion of similar spectral type and 4-8 L Q can 
only be a red (sub-)giant of slightly lower initial mass than 
Arcturus. Note that a higher mass is not possible because 
of evolutionary reasons: in such case, the companion would 
have been much brighter than the primary, under the as- 
sumption that they are coeval. 

A quantitative estimate of the allowed mass ra- 
tio and spectral types can be made using synthet- 
ical isochrones/evolutionary tracks. Using the Padova 
datab ase of stellar evolutionary tracks <|Girardi et alJ 
l2000ft fo r z = 0.004, and the luminosity, T e g and log g 
of lDecinl l|2000j) . we find for the primary (and the possible 
system) an age range from about 10 985 to 10 10 - 25 years 
(7 to 17 billion years, actually limited by the age of the 
universe). In Fig.^J we present luminosity as a function 
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Fig. 14. Luminosity as a function of initial mass for a 
system of 10 9 85 and 10 10 25 years old. For the companion, 
we indicate the possible range in luminosity The two peaks 
in luminosity on the right hand side correspond to the 
RGB and AGB phase respectively. 



of initial mass for the lower and upper limit on the age of 
the system. 



Table 7. Primary and companion stellar parameters for 
age-limits 2 . 



Log(Age) 




9.85 


10.25 


Component 


A 


B 


A 


B 


L (L Q ) 


195 


3.8-8 


195 


3.8-8 


T cff (K) 


4574 


5050-6300 


4450 


4950-5050 


log g 


1.69 


3.5-4 


1.52 


3.5 


Minitial 


1.03 


0.99-1.02 


0.805 


0.799-0.802 


M ac tual 


1.02 


0.99-1.01 


0.796 


0.796 


Spect. Type 


K0III 


G5IV-F7IV 


K1.5III 


G5IV 
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In Table [7| we present the stellar parameters for pri- 
mary and companion at the two edges of the age interval 3 . 
From this ta ble, it i s clear that the stellar parameters de- 
termined bv iDecinl l)2000(l correspond to a fairly old sys- 
tem, but the values found for the 7-billion-year system are 
still compatible with the literature on a Bootis. Note that 
the lower mass in the "old" solution corresponds to the 
mass estimate implied by the measured radius and the 
derived gravity, which is close to 0.8 M©. Furthermore, 
the V— K colour (1.98) one obtains for the companion 
when combining our K-band solution with the Hipparcos 
V-ba nd solution corresp onds to a spectral type of G4IV 
(e.g. iBessell et al.lll998l) . exactly what can be expected 
from Table [7] and our preference for a high-age solution! 
To obtain a K-band contrast of 50, the companion would 
need to be about 2.7 mas large. This value is entirely com- 
patible with the evolutionary model of the companion put 
at a distance of 11.26±0.09pc (derived from the parallax) 
and with the diameter com puted with the empiri cal sur- 
face brightness relation from lKervella et all l|2004lh which 
yields 2.72 ± 0.04. It is compatible as well with the bi- 
nary parameters reported in Table H3 as recomputing the 
binary model with 02 = 2.7 mas does not lead to a no- 
ticeable change in the solution. 

In the lower panel of Fig. we present the abso- 
lute flux of this hypothetical companion together with the 
residuals we obtain by subtracting a model for the pri- 
mary from the observed photometry 4 . Clearly, the com- 
panion would not be detectable in these residuals, except 
maybe for the B and V bands (U should be treated with 
much care because it is highly surface-gravity sensitive), 
which show significant excess, compatible with the sug- 
gested 2. 7- mas G4IV companion. 

Note that we found the current secondary mass to be 
almost equal to the primary mass. A mass ratio so close to 
unity, unlikely as it seems, is compatible with statistical 
binary mass ratio studies showing a bimo dal distributio n 
with one peak toward a ratio of 1 (e.g. lTrimblelll97^ . 
The difference in luminosity is the result of the companion 
being slightly behind in evolution. 

5.2.5. The orbit 

Combining all of these results, two plausible types of orbit 
remain: (1) a very narrow system (separation of just a few 
AU) seen face-on or (2) a very wide system (period of the 
order of centuries, or larger) seen nearly edge-on. 

The first solution does require us to discard the 
PA/relative motion detected in the Hipparcos data, 

3 The actual T e g determined bv lDecinl l)200fl) is beyond the 
grid-limit, but by extrapolation one can see that this probably 
does not influence our results. 

4 The diameter of the primary was chosen in such a way 
as to obtain positive residuals (except for the J-band which is 
notoriously sensitive to the water column density in the earth's 
atmosphere). This required a slightly smaller diameter than 
derived in Sect. l4.l"!?l i.e. 20.0 mas. 



which is acceptable g iven the doubts expressed by 
ISoderhielm fc Mignardl (^998) in the re-analysis of the 
data. It would however be compatible with the periods 
found in radial velocities and astrometry of the order of a 
few years. To explain the low observed radial velocity vari- 
ations, the inclination must be very close to face-on, or the 
companion must have a mass of just a few Jupiter masses, 
which is not compatible with the observed luminosity. 

The second solution would require a nearly edge-on 
disk to produce the small projected separation, and the 
long orbital period would explain the non-detection of 
high-amplitude radial velocity and astrometric variabil- 
ity. Clearly, it is not possible to combine all results into 
one consistent explanation. 

6. Conclusions and Outlook 

We have presented a new set of narrow-band near-IR inter- 
ferometric observations of the K2 giant Arcturus, obtained 
to test the applicability of K giants and K giant models to 
the calibration of high-accuracy infrared interferometric 
observations. A comparison with state-of-the-art stellar 
atmospheres failed to explain the data: the residuals show 
clear systematics. These are independent of calibration 
and none of our hypotheses about the stellar atmosphere 
(extension, spots) can explain the data. A binary model, 
with a sub-giant as companion, does provide a good fit to 
the data. A thorough literature study reveals that there 
is ample, but inconclusive and inconsistent evidence for a 
companion. 

6.1. Implications on calibration 

If the calibrator one uses is an unknown binary, or the 
effects reported in this paper are due to some other un- 
known characteristics of K giant stars, then the error on 
the visibility could be increased by a factor as large as 
v/Xr = \/5.6 Ri 2.5, if a standard model atmosphere is 
used to represent the calibrator. For single-mode inter- 
ferometers such as FLUOR or VINCI, it means that the 
standard error on the visbility would jump from 1-2 % to 
2.5-5%. Further investigation is urgently required to set- 
tle this major issue. This should include both new observa- 
tions of other "normal" K giants to check the uniqueness 
of the problems reported here, and the re-observation of 
Arcturus to pin-point with more certitude the nature of 
this star. 

6.2. Prospects for Arcturus 

Unfortunately, our dataset which was aimed at deter- 
mining the wavelength dependence of the diameter of 
Arcturus (A), does not allow a full characterization of 
the (still hypothetical) system. New attempts at imaging 
with the most recent AO systems might resolve the sys- 
tem, though the brightness of the primary will seriously 
complicate the observations. Another possibility would be 
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a new interferometric dataset with two telescopes on sim- 
ilar or shorter baselines, but obtained over a wide range 
of azimuth angles, such that the companion would pass 
through several of the fringes, providing a much clearer 
signal than was possible with the present essentially snap- 
shot data set. Yet another possibility would be an inter- 
ferometric data set with three telescopes (e.g. the present 
IOTA), using two short baselines to lock on the fringes, 
and one long baseline on which the large-diameter primary 
would have weak fringes but the small-diameter secondary 
would have relatively strong fringes, again observing over a 
large range of position angles so that the secondary passes 
through many fringes, for a clear signature. 
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